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ABSTRACT 

We report on the long term stability of the millisecond oscillations observed with 
the Rossi X-ray Timing Explorer (RXTE) during thermonuclear X-ray bursts from the 
low mass X-ray binaries (LMXB) 4U 1728-34 and 4U 1636-53. We show that bursts 
from 4U 1728-34 spanning more than 1.5 years have observed asymptotic oscillation 
periods which are within 0.2 /isec of each other, well within the magnitude which 
could be produced by the orbital motion of the neutron star in a canonical LMXB. 
This stability implies a timescale to change the oscillation period of > 23, 000 years 
in this system, and suggests a highly stable process, such as stellar rotation, as the 
mechanism producing the oscillations. For 4U1636-53, which has an orbital period of 
3.8 hours, we show that offsets in the asymptotic oscillation periods from three different 
bursts can be consistently interpreted as due to orbital velocity of the neutron star with 
vsmi/c ~ 4.25 x 10~ 4 . An updated optical ephemeris for the epoch of maximum light 
from V801 Arae would provide a strong test of this interpretation. We discuss the 
constraints on the X-ray mass function which can in principle be derived using this 
technique. 

Subject headings: X-rays: bursts - stars: individual (4U 1636-53, 4U1728-34) stars: 
neutron - stars: rotation 



1. Introduction 

Millisecond oscillations in the X-ray bright- 
ness during thermonuclear bursts have now been 
observed from 6 low mass X-ray binaries (LMXB) 
with the Rossi X-ray Timing Explorer (RXTE) 
(see Strohmayer et al. 1996; Smith, Morgan & 
Bradt 1997; Strohmayer et al. 1997; Zhang et 
al. 1996; Zhang et al. 1997; and Strohmayer, 
Zhang & Swank 1997). The presence of large am- 
plitude oscillations near burst onset, combined 
with spectral evidence for localized thermonu- 
clear burning suggests that the oscillations are 
caused by rotational modulation of thermonu- 
clear inhomogeneities (see Strohmayer, Zhang & 
Swank 1997). 

The accretion-induced rate of change of the 
neutron star spin frequency in a LMXB is ap- 
proximately 
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where rhn, M x , r acc and Z45 are the mass accre- 
tion rate in units of 10 17 g s _1 , the neutron star 
mass in solar units, the characteristic radius of 
the inner accretion disk in km, and the stellar 
moment of inertia in units of 10 45 g cm 2 , respec- 
tively. If the millisecond oscillations observed 
in the X-ray brightness from thermonuclear X- 
ray bursts with the Rossi X-ray Timing Explorer 
(RXTE) are produced by rotational modulation 
of the burst flux, then the Doppler corrected 
frequencies should be stable at the better than 
Au = 0.001 Hz level over a hundred years or so. 
The Doppler shift due to orbital motion of the 
binary can produce a frequency shift of magni- 
tude 
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where M x , M c , Ph r , v and i are the neutron 
star mass, the companion mass (both in solar 
units), the orbital period in hours, the magnitude 
of the neutron star orbital velocity, and the sys- 



tem inclination angle, respectively. For canonical 
LMXB system parameters this doppler shift eas- 
ily dominates over any possible accretion-induced 
spin change on orbital to several year timescales. 
Thus, the level of observed stability in oscillation 
periods from burst to burst provides a method to 
further test the rotational modulation hypoth- 
esis. For example, if oscillation period shifts 
larger than can plausibly be produced via or- 
bital motion are observed this would tend to cast 
doubt on the spin modulation interpretation. On 
the other hand, if the burst oscillation frequen- 
cies remain stable over long timescales, reveal- 
ing a signature of binary motion, then it will 
both support the rotational interpretation and 
become possible to use the observed frequency 
shifts to constrain the neutron star binary mass 
function in systems which have observed burst 
oscillations. 

In this Letter we investigate the long term pe- 
riod stability of burst oscillations in two LMXB 
sources, 4U 1728-34 and 4U 1636-53. We show 
that in 4U 1728-34 bursts separated in time by 
about 1.6 years have shortest observed oscilla- 
tion periods, what we refer to as asymptotic pe- 
riods, within 0.2 /xsec of each other, well within 
the range of shifts which could result from the 
system's binary motion. For 4U 1636-53, which 
has a known orbital period of 3.8 hours, our 
time baseline is shorter, however, the orbital pe- 
riod allows us to show that the oscillation pe- 
riod shifts observed from three different bursts 
can be consistently interpreted in terms of those 
produced by binary orbital motion with reason- 
able values for the component masses and sys- 
tem inclination. If the relative orbital phases 
when the bursts ocurred can be converted to ab- 
solute phases, for example, with an updated op- 
tical ephemeris of the epoch of maximum light, 
then this would provide a test of the doppler 
shift interpretation, and if confirmed would en- 
able constraints on vsini/c and thus the X-ray 
mass function to be derived. 
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2. Long Term Frequency Stability in 4U1728- 
34 

We had observations of 4U1728-34 with RXTE 
in February, 1996 and again in September, 1997 
(see Strohmayer et al. 1996 for a summary of the 
February, 1996 observations). Using data from 
these two observations we can compare the oscil- 
lation periods during bursts over a span of about 
1.6 years. For all the burst data reported here we 
had either 125 /j,sec (1/8192 s) resolution binned 
data or event mode data with the same temporal 
resolution. 

The oscillations at 2.75 ms (363 Hz) observed 
during bursts from 4U 1728-34 are not strictly 
coherent (see Strohmayer et al. 1996). In some 
bursts the period is observed to evolve from 
a high of about 2.762 ms near burst onset to 
about 2.747 ms during burst decay. Due to 
the episodic nature of the oscillations not all 
bursts show detectable oscillations over this en- 
tire range. Strohmayer et al. (1997) have argued 
that this frequency evolution is caused by the in- 
crease in the scale height of the thermonuclear 
burning layer on the neutron star surface and 
subsequent conservation of angular momentum 
of the thermonuclear shell. In many bursts which 
show oscillations the oscillation period appears 
to reach a nearly coherent, asymptotic limit as 
the burst decays away. In the context of the spin 
modulation hypothesis this limit represents the 
actual spin period of the bulk of the neutron star. 

We selected for detailed comparison a pair 
of bursts from the February, 1996 observations 
(bursts 4 and 5, respectively, from Strohmayer, 
Zhang & Swank 1997) and one from the Septem- 
ber, 1997 data. Here we refer to these bursts as 
bursts 1, 2 and 3, in time order. We selected 
these bursts because they showed significant os- 
cillations over the longest time intervals during 
the bursts and the oscillation period during the 
burst decay reached a stable, coherent limit. In 
Figure 1 we compare the dynamic power spec- 
tra of the bursts detected on Feb. 16, 1996 at 



06:51:07 UTC and Sep. 22, 1997 at 06:42:51 UTC 
(bursts 2 and 3 in Table 1). We only show two 
of the three bursts since the pair of bursts ob- 
served on Feb. 16, 1996 were nearly identical 
in their oscillation properties. The figure shows 
contours of constant power spectral amplitude 
and they have been shifted in time for clarity. 
The dynamic power spectra were computed from 
2 s intervals with a new interval starting every 
1/8 s. The leftmost contours are for the Feb., 
1996 burst. The frequency evolution from low 
to high is clearly evident in both bursts, and the 
range of observed frequencies and the highest ob- 
served frequency are very similar. Note that the 
oscillation frequencies in both bursts reach a sta- 
ble upper limit, what we will call the asymptotic 
frequency or period. 

Since the rotational modulation hypothesis 
suggests that the shortest observed period is the 
underlying stellar spin period we performed an 
epoch folding period search analysis using only 
the portions of all three bursts after which the 
frequency has stabilized to see how closely these 
asymptotic frequencies agree. Figure 2 shows 
the resulting \ 2 plots from the folding analysis 
function of barycentric period. To estimate 
the oscillation periods and uncertainties from the 
epoch folding we computed the centroids P cen 
and standard deviations ap of each x 2 peak using 
the relations 

p _ 2~Zi XjPi 2 _ XjjPi ~ Pcen) 2 /„\ 

r- cen - a P — 2 , \o) 

Xi 2—ii Xi 

where i runs over each bursts x 2 peak from the 
epoch folding analysis. Table 1 summarizes the 
derived asymptotic periods and uncertainties for 
each burst. As can be seen from the inferred pe- 
riods and uncertainties there is no significant evi- 
dence that the observed asymptotic periods from 
the three bursts are different. Using the mea- 
sured centroids as the best estimator of the peri- 
ods for each burst, the implied period difference 
over the 1.6 yr timespan is about 0.19 fisec. In 
terms of a timescale to change the period this cor- 
responds to r > P/P = PAT/AP = 2.3 x 10 4 yr, 
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and implies a limit on any orbital doppler shift 
AP/P = v or bSva.i/c < 6.9 x 1CT 5 , well within the 
shift which could be produced by orbital motion 
of the neutron star in a typical LMXB. 

3. Burst Oscillation Frequencies in 4U 
1636-53 

X-ray brightness oscillations during X-ray bursts 
at 1.72 ms (581 Hz) were discovered in 4U 1636- 
53 by Zhang et al. (1996). The 3.8 hr orbital 
period of 4U 1636-53 is known from the observed 
optical periodicity of the optical companion V801 
Arae (see van Paradijs et al. 1990; Smale & 
Mukai 1988; and Pedersen, van Paradijs & Lewin 
1981). Since the orbital period is known the rel- 
ative phases at which bursts occurred can be de- 
termined. One can then compare the observed 
oscillation periods from different bursts and de- 
termine if any observed changes can be consistent 
with an orbitally induced doppler shift. In par- 
ticular, if at least three bursts are available with 
measured oscillation periods then one can try to 
solve the following set of equations; 

P t . =P - APcos(<&, + fa) ■ (4) 

Here, Pt j and fa. are the observed asymptotic os- 
cillation periods and relative orbital phases, re- 
spectively, for bursts which occurred at tj, and 
Pq, AP, and 4>o are the barycentric oscillation 
period when the neutron star transits the line 
of sight, the magnitude of the doppler induced 
period change, and an initial phase offset, re- 
spectively. With at least three different oscilla- 
tion period measurements during bursts it may 
be possible to determine a set of values for the 
three parameters Pq, AP, and (f)Q which are con- 
sistent with binary motion. The orbital veloc- 
ity v and system inclination i are related by 
AP/Pq = vsmi/c. 

For 4U 1636-53 we now have 3 different bursts 
spanning a time interval of slightly more than 
a day. To determine the asymptotic oscillation 
periods for each burst we performed a similar 



epoch folding analysis on these bursts as those 
from 4U 1728-34 described above. Table 1 sum- 
marizes information on the occurrence times, rel- 
ative orbital phases (from burst 1) measured at 
the solar system barycenter of 4U1636-53 at the 
time of occurrence, and the barycentric asymp- 
totic period observed during the decaying portion 
of each burst. Figure 3 displays the resulting \ 2 
plots from the epoch folding analysis for each of 
the three bursts. 

Bursts 1 and 3 occurred approximately half 
an orbit apart from each other and note that 
these bursts had asymptotic periods which are, 
within the uncertainties, consistent with each 
other. Burst 2 occurred roughly midway in or- 
bital phase between bursts 1 and 3 but had a 
significantly shorter period by about 0.74 ^sec. 
Thus, a plausible scenario is that burst 1 oc- 
curred near the time of superior conjuction of 
the neutron star (inferior conjunction of the sec- 
ondary). At this phase, as well as half an orbit 
away, the neutron star is near transit and the or- 
bital component of its velocity along the line of 
sight is nearly zero. This can explain the consis- 
tent periods measured during bursts 1 and 3. Fi- 
nally, burst 2 occurred 0.327 in phase from burst 
1 and thus with a significant fraction of its total 
orbital velocity along the line of sight to produce 
the strong blue shift to shorter period. 

This scenario gives vsmi/c « 4.25 x 10~ 4 to 
account for the observed period shift between 
burst 2 and bursts 1 and 3. Current understand- 
ing of the mechansim underlying the optical mod- 
ulations at the orbital period suggests that the 
epoch of maximum light be identified with su- 
perior conjunction of the optical secondary (see 
Smale & Mukai 1988). Van Paradijs et al. (1990) 
give an ephemeris for the epoch of maximum 
light for V801 Arae based on a compilation of 
observations from as early as July, 1980 to as 
late as May, 1988. However, the uncertainties in 
projecting this ephemeris forward to the epoch 
of burst 1 are such that it does not provide a sig- 
nificant test of the orbital doppler interpretation, 
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but, with an updated optical ephemeris a much 
more restrictive test can be made. 

Since the orbital motion hypothesis cannot be 
rejected with the present data we show in fig- 
ure 4 the contours of constant sini in the mass 
plane for the appropriate orbital period of 3.8 hr 
and the vsini/c suggested by the period shifts in 
the bursts. For the 3.8 hour orbital period and a 
main sequence companion, empirical mass-radius 
relations would suggest a mass of about 0.36 
M sun for V801 Ara (see Patterson 1984; Smale 
& Mukai 1988). This is a bit lower in mass than 
that suggested by the inferred v sin i/c from the 
bursts (see Figure 4) , but given the uncertainties 
in the period measurements as well as the uncer- 
tain effects of X-ray heating on the secondary, it 
is still within a reasonable range for the observed 
shifts to be plausibly produced by orbital mo- 
tion. Both the observation of additional bursts 
as well as an updated optical ephemeris for V801 
Ara can provide a careful test of the orbital hy- 
pothesis for the observed burst oscillation period 
shifts. 

4. Discussion 

The episodic nature of the oscillations dur- 
ing bursts does introduce some uncertainty into 
what is the "highest" observed frequency dur- 
ing a burst. It is possible that in some bursts 
the oscillations are not strong enough to be de- 
tected at late times in the burst and therefore 
the highest frequency may not be observed in all 
bursts. Partly this can be mitigated by compar- 
ing bursts which show similar overall oscillation 
properties, as we have endeavored to do here, 
but in order to fully overcome this one simply 
needs the weight of evidence from a larger sam- 
ple of bursts. In particular, with a large enough 
sample to cover most of the orbital phase space 
the signature of orbital doppler shifts should be- 
come fairly transparent, or not, since the mag- 
nitude of the frequency offsets should be limited 
by the magnitude of the binary orbital velocity, 
and an approximately equal number of redshifts 



and blueshifts should be observed. 

The long term stability of the highest millisec- 
ond oscillation frequencies observed in thermonu- 
clear bursts from 4U1728-34 and 4U1636-53 pro- 
vides a strong argument in favor of a higly sta- 
ble clock, such as stellar rotation, setting the ob- 
served oscillation frequency. Regardless of the 
mechanism, any oscillation period will suffer or- 
bital doppler effects. The limits on the period 
offsets from bursts spanning 1.6 years in 4U 1728- 
34 indicates that the intrinsic period which sets 
the asymptotic period during bursts can change 
on a timescale no shorter than r = PAT/ AP ss 
2.3 x 10 4 yrs. This timescale is longer than sim- 
ilar timescales for many known X-ray pulsars, 
and is also longer than the characteristic time 
to change the thermal state of the neutron star 
surface ocean (see Bildsten, et al. 1998). Thus if 
oscillation modes sensitive to the thermal state, 
such as g- modes, were the cause of the oscilla- 
tions, they would not be expected to be stable 
over such a long timescale. If analysis of addi- 
tional bursts continues to support this interpre- 
tation, then it will become possible to use os- 
cillation periods observed in different bursts to 
place constraints on the masses of the compo- 
nents in LMXB, thus long pointed observations 
that collect many bursts are well justified given 
that they could lead to a determination of the 
mass function for a larger sample of systems. In 
addition, constraints on v sin i jc derived from dif- 
ferent bursts provides a method to conduct more 
sensitive searches for the millisecond X-ray pul- 
sar in the persistent, accretion-driven flux, which 
should be present at some level in most LMXB. 

We thank Nick White, Glenn Allen and Mike 
Stark for helpful discussions and comments on 
the manuscript. 
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Table 1: Asymptotic Oscillation Periods in bursts from 4U1728-34 and 4U1636-53 



Source 


Tburst (UTC) 


P cen (ms) 


op (/j,sec) 


4>rel(TBD) 


4U 1728-34 


burst 1 


2/16/96 at 06:51:13 


2.74774 


0.40 


NA 


burst 2 


2/16/96 at 10:00:49 


2.74757 


0.39 


NA 


burst 3 


9/22/97 at 06:42:56 


2.74755 


0.39 


NA 


4U 1636-53 


burst 1 


12/28/96 at 22:39:24 


1.71970 


0.43 


0.0 


burst 2 


12/28/96 at 23:54:03 


1.71896 


0.39 


0.327 


burst 3 


12/29/96 at 23:26:47 


1.71981 


0.38 


0.517 
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Fig. 1. — Dynamic power spectra computed from two different bursts from 4U 1728-34 separated in time 
by 1.6 yr. Shown are contours of constant power spectral density. The contours have been offset from 
each other for clarity. Note that the range in frequency of the oscillations as well as the highest observed 
frequency are very similar. The burst from 2/16/96 at 10:00:49 UTC is on the left, that from 9/22/97 at 
06:42:56 UTC in on the right. 
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Fig. 2. — Results from the x 2 epoch folding analysis for the three bursts from 4U 1728-34. The bursts 
are arranged in time order from bottom to top (burst 1 at bottom to burst 3 at top). See Table 1 for the 
measured period centroids and uncertainties. 
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Fig. 3. — Results from the x 2 epoch folding analysis for the three bursts from 4U 1636-53. The bursts 
are arranged in time order from bottom to top (burst 1 at bottom to burst 3 at top). See Table 1 for the 
measured period centroids and uncertainties. 
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Fig. 4. — Contours of constant sini for a binary system with orbital period 3.8 hr and vsmi/c = 
4.25 x 10 -4 as suggested by the observed period offsets in bursts from 4U1636-53. This result should not 
yet be taken as a constraint on the system masses in 4U 1636-53, rather it only suggests that the orbital 
motion is a plausible explanation for the observed period shifts. 
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